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Abstract 

We study a three dimensional Z(3)-symmetric effective theory of high temperature QCD. 
The exact lattice-continuum relations, needed in order to perform lattice simulations with 
physical parameters, are computed to order ©(a^) in lattice perturbation theory. Lattice 
simulations are performed to determine the phase structure of a subset of the parameter 
space. 
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1 Introduction 



At high temperature, QCD matter undergoes a deconfinement transition, where ordinary had- 
ronic matter transforms into strongly interacting quark-gluon plasma [H EJ [3] . In the absence 
of quarks, Nj = 0, the transition is a symmetry-breaking first order transition, where the order 
parameter is the thermal Wilson line [U [5] . The non-zero expectation value of the Wilson line 
signals the breaking of the Z(3) center symmetry of quarkless QCD at high temperatures. 

The transition has been studied extensively using lattice simulations [El O |8| . Thermody- 
namical quantities, condensates and various correlators can be measured on the lattice and the 
equation of state can be estimated. This approach, however, becomes computationally exceed- 
ingly expensive at high temperatures, and thus cannot be applied to temperatures T above 
~ 5Tc. The complementary approach has been to construct perturbatively effective theories, 
such as electrostatic QCD or EQCD, using the method of dimensional reduction [9l [lOl E] to 
quantitatively describe high temperature regime of QCD [T21 ESI EH ES]- In the dimensional 
reduction procedure, however, one expands the temporal gauge fields around one of the Z(3) 
vacua and thus explicitly violates the center symmetry. The range of validity of these theories 
therefore ends for T below ~ 5Tc, where the fluctuations between different vacua become im- 
portant. There have also been several attempts to the build models for Wilson line, respecting 
the center symmetry [l6l E3 EH EHl EOl EI] • These models give a qualitative handle on the 
transition but cannot be perturbatively connected to QCD. 

As a unification of these strategies, an effective field theory of high temperature QCD respect- 
ing the Z(3) center symmetry has been constructed in [22]. At high temperatures, the effective 
theory reduces to EQCD guaranteeing the correct behavior there, but the model still preserves 
the center symmetry. The effective theory is further connected to full QCD by matching the 
domain wall profile separating two different Z(3) minima. 

Being a three dimensional model, the new theory relies on the scale separation between the 
inverse correlation length and the lowest non-zero Matsubara mode, which is still modest at Tc 
[231 124j . Thus, one hopes that the range of validity of this theory would extend down to Tc- 
The effective theory is a confining one, so perturbative analysis breaks down. Non-perturbative 
methods, i.e. lattice simulations, are thus needed to find out the physical properties, such as 
correlation lengths, condensates, and most importantly the phase structure of the theory, to test 
its regime of validity. 

The effective theory is super-renormalizable, and thus the connection between the continuum 
MS and lattice regulated theories can be obtained exactly to the desired order in the lattice 
spacing a. The matching of the parameters of the Lagrangian to order 0{a^), which is needed 
in order to perform simulations with MS scheme parameters and to obtain physical results, 
requires a two-loop lattice perturbation theory calculation. The one-loop terms remove any 
linear 1/a divergences, while two-loop terms remove the logarithmic log(l/a) divergences and 
the constant differences in the mass terms of the theory. In addition, the condensates have lattice 
spacing dependence and constant differences between the two schemes and can be calculated to 
order 0{a^) by performing a two-loop calculation for operators up to cubic order and a four- 
loop calculation for the quartic condensates. In this paper we perform the needed two-loop 
calculations. 

This paper is organized as follows. In Sections 2 and 3 we define the theory in continuum 
MS regularization and on the lattice, respectively. In Section 4 we study the phase diagram of 
a subset of the parameter space of the theory. Details of the matching between the continuum 
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and lattice theories are given in the appendices. 



2 Theory 



The theory we are studying is defined by a three dimensional continuum action, which we 
renormalize in the MS scheme 



S = I d3-2^x I^Tr^ + Tr (a^^A^) + Vo{Z) + Vi{Z)^ 



where 



diAj 



djAi + ig3[Ai,Aj] 



di - igslAi 



(1) 



(2) 
(3) 



and Z is a 3 X 3 complex matrix, which in the limit e — > has dimension dimZ = V GeV. 
The gauge fields Ai are Hermitean traceless 3x3 matrices and can be expressed using gener- 
ators of SU(3), Ai = AfT", with TiT^T^ = ^6'^^ . The covariant derivative is in the adjoint 
representation. The potentials Vq, the "hard" potential, and Vi, the "soft" potential, ar^: 



Vo{Z) = ciTr[Ztz] + 2c2Re(Det[Z]) + C3Tj:[{Z^ Z)% 
Vi{Z) = diTr[M^M] + 2(i2Re(TV[M3]) + d3Tr[(M^M)2] 



(4) 
(5) 



where M = Z — |Tr[Z]l is the traceless part of Z. Here, the gauge coupling has a positive 
mass dimension diml^g] =GeV, making the theory super-renormalizible. Because of the super- 
renormalizibility, the coefficients 02,03,^25 and ds are renormalization scale independent and 
only the mass terms ci and di acquire a scale dependence in the MS renormalization scheme. 
The scale dependence in the mass terms arises from a two-loop calculation and has the form: 



ci(m) 
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(6) 
(7) 



where A is a constant specifying the theory and \i is the MS scale parameter. The coefficients 
Cj, dj, and gj, are matched to the parameters of full thermal QCD by imposing the condition 
that the theory reduces to EQCD at the high temperature limit, and that the theory reproduces 
the domain wall profile of full QCD [22]. This defines a subset of parameter values (with a 
limited accuracy due to perturbative matching), for which the theory describes thermal qceI. 
However, in this paper we consider the model in general, and do not restrict ourselves only to 
the physical region. 

The action is defined only for the number of colors A'^c = 3, but for generality, we give some of 
the perturbative results for any Ac. For analytic calculations the scalar field Z can be expanded 
around the vacuum: 

Z = \\\(<^^iX)^^iM^iA), (8) 



^Our notation is obtained from that in [22] by scaling with gs: Ai QsAi, Z g:iZ, ci ci, C2 — > 33 ^C2, 
C3 93^ C3, ci di, C2 ^ 93 ^d2, and C3 93^ ds 

•^In the matching, the hard potential is parametrically larger than the soft potential, explaining the terminology. 
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where (f> and % are real scalars and H and A are Hermitean traceless matrices. Fields H and A 
can be written with the generators of the SU(3) group, H = H^T^ and A = ^"T", where iJ" 
and A"" are real scalars. 

The action is invariant under local gauge transformations, with Z transforming in the adjoint 
representation: 

Ai{x) G{x) (^Ai{x) - G-\x), (9) 

Z{x) — > G{x)Z{x)G-\x), (10) 

where G{x) gSU(3). In addition to this, there are further global symmetries in the potentials. 
The potential Vq is invariant in global SU(3)xSU(3) transformations 

Z{x) — ^ LZ{x)R, (11) 

where L and R are SU(3) matrices. The potential Vi is invariant under Z(3) transformations 

M zM, (12) 

where z = e*^'^"/'^, which generalizes into a U(l) symmetry if d2 = 0. This implies that in the 
presence of the both potentials (with non-zero coefficients), the overall global symmetry of the 
Lagrangian is Z ^ zZ. 



3 Lattice action 



In order to perform non-perturbative simulations, the theory has to be formulated on the lattice. 

On the lattice, the scalar field Z lives on the sites of the lattice, and the gauge fields Ai are traded 
for link variables Ui, which arc elements of SU(A''c) and live on the links connecting adjacent 
sites. The lattice action corresponding to the continuum theory can be written as ^ = Sw + Sz, 
where 



1 



1 _ —ReTriU, 



/if J 



x,i<j 

is the standard the Wilson action with the lattice coupling constant 

2Nr, 



(13) 



2 ' 



corresponding to a lattice spacing a. The continuum limit is taken by /? — ^ oo, and there the 
Wilson action reduces to the ordinary pure gauge action. 

The kinetic term, Tr (L>jZ^^Z)jZ), is discretized by replacing the covariant derivatives by 
covariant lattice differences. Then the scalar sector of the action reads: 

2Nr 



Sz = 2[-^) ^ReTr Z^ Z - Z\x)Ui{x)Z{x + i)t//(a 



+ 



+ 



/ ^{ciTi[Z^Z]+2c2ReDeiZ + c^lh:[{Z^Zf\^ 



(14) 



3 



where Cj, di are dimensionless numbers, and M and Z are dimensionless Nc x Nc complex matri- 
ces. Only the mass terms ci and di require non-trivial renormalization and all the other terms 
can be matched to order 0{a^) on tree-level by simply scaling with g^: 



Z 

C3 = 



ales, 



M = gsM 
d2 = gld2 
ds = gids- 



(15) 
(16) 
(17) 



For the mass terms, renormalization has to be carried out, so that the physical masses of the 
fields are the same in both regularization schemes. A two-loop calculation gives (the details of 
the calculation and the definitions of the numerical constants are given in the appendix): 
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and 
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^ci - 644 + -^4c3 + -dl + I J [log/3 + 0.08849] } + 0(/3-i). (19) 

Here, we have set the renormalization scale to be /i = g| in Eqs. ^ and ([7]), and denote 
Cl = c\{g1) and di = di{g1). By making this choice, we get the logarithmic term to be a 
function of the lattice coupling constant (3. There are also higher order corrections (corrections 
of order C'(/3~^) corresponding to order 0{a) in lattice spacing), but their effect vanishes in the 
continuum limit. 

Various operators also need to be renormalized on the lattice in order to convert their expec- 
tation values to continuum regularization. The Zs-symmetry protects the lowest dimensional 
condensate {g^ TvZ) from acquiring any additive renormalization, while a two-loop calculation 



4 



gives for the quadratic condensates: 



{gfT,Z^Z)y^={i:,Z^Z)a 



2iVc(iV2-l) log/3 + C+ — -5 



(53-'TVMtM)j^ = (TrMtM), 



{TvZ^Z)a - [0.3791/3 + 0.3040(log /? + 0.66796)] + 0{(5-^] 
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(TrZtTrZ),-|-|-/3 + 0(r'), 



(20) 

(21) 
(22) 



{g^'TrZ^TiZ)^ 

where the subscript a denotes the lattice regularization. For the cubic condensates we get: 
(2ReDeti')a 

' (iog/3 + c) + o(r' 



?3-32ReDetZ)Hs 
1 
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The effect of subtraction of the divergences can be seen in Fig.l. The renormahzation of the 
quartic operators to order 0{I3^) would require a four-loop calculation, which we do not perform 
here since they are not measured at this stage. 
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Figure 1: The effect of subtraction of the divergences in (TrZ^Z) and (TrA/^A/) in a fixed physical 
volume with di = 6.6 and — 0.01. On the left panel, the elfect of subtracting the divergent parts 
of (TrZ^Z) is plotted. The constants ki, K2 and K2 are the coefficients of the linear, logarithmic, and 
constant differences between lattice and MS regularizations form equation (|^D|) . On the right panel: 
the continuum limit of the condensates. Notice the negative values of the quadratic condensates in the 
symmetric phase. 
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4 Phase diagram of the soft potential 



A simpler model is obtained from the original theory by setting Cj = 0. In this model, the trace 
of Z decouples and can be integrated over as a free scalar field. The relevant degree of freedom 
is thus a traceless complex matrix M, or two traceless hermitian matrices H. and A. This can 
be viewed as a natural generalizatiorcl of EQCD to complex values of the adjoint higgs field A!^. 
The simpler model is defined by the action: 

5- = y d^x ^Tri^^j. _^TrAM^AAf + diTrMtM + 2d2Re(Tr[M3]) + d3Tr(M^M)2 (25) 

= y d^xj^Tri^^ _^ TrDiADiA + TrDiHDiH + diTiA'^ + diTri?^ 

+ 2d2Tr[if3 - SHA^] + d3Tv[H^ + A^ + AH'^A^ - 2HAHA]Y (26) 

If the cubic term ^2 is zero, the Lagrangian is invariant under a U(l) global symmetry 
M gM, g gU(1). The breaking of the symmetry is signalled by a local order parameter: 



A = ^(1^^3)2 + (T¥//3)2. (27) 

This operator remains a valid order parameter after the renormalization since it has no additive 
renormalization, if ^2 = 0. In the symmetric phase A is strictly zero and in the broken phase the 
order parameter obtains a non-zero vacuum expectation value, while the two phases are separated 
by a first order transition. In the broken phase (TrM^M) is larger than in the symmetric phase. 
After the inclusion of the cubic term, A is no longer strictly an order parameter, since the U(l) 
symmetry is explicitly broken. However, the first order transition remains and is accompanied 
with a significant discontinuity in A and (TrM^M). 



4.1 Perturbation theory 

In the limit of small d^/g^, the transition becomes very strong, and we expect a semiclassical 
approximation to produce the correct behavior of the critical line [25\ I26j . We parametrize a 
constant diagonal hermitian background field in a fixed Landau gauge as follows: 

/ q+P \ 
(M) = 2pT3 + 2V3qTs = q-p , (28) 

\ -2q J 

where p and q are real scalars with dimensions of g^. 

Lattice simulations suggest that the A — > —A symmetry is not broken spontaneously at any 
non-zero value of d2, so that it is sufficient to consider only hermitian background fields. Using 
this parametrization, the 1-loop effective potential Vi{di,d2,d3;p,q) can be calculated: 

EQCD with gauge group SU(3), there is only one linearly independent quartic gauge invariant opera- 
tor namely Tr^lg. In the complex case, however, there are four different Za-symmetric operators: Tr(M^M)^, 
(TrMtM)^ TrfM^M^MM] and Tr[Mt^]Tr[M^]. In the case of unitary M, i.e. in the minimum of the hard 
potential, these operators collapse into a single one. However, since there is no such restriction in our model, the 
operators are linearly independent. From these operators, we choose to include only the one appearing in the 
original theory, Tr(M'l"M)^. 
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Vi{di,d2, ds; q,p) = (2p' + 6g2)di + 36q{p^ - q^)d2 + -{2p' + 6^2)2^3 



-^(8H' + |p-3g|^ + |p + 3g|3)5f 
^{2 [di + 3(p - q)d2 + 2(p2 + 3g2)d3 
+2 [di + 6gd2 + 2(/ + 3g2)t^3 
+2 [di - 3(p + g)d2 + 2(p2 + 3g2)d3 



3/2 



i3/2 



i3/2 



+ 



+ 



di + 4(p2 + 3g2)d3 _ 2 J3(p2 + 3g2)(i2 + 18g(p2 _ ^^2)^2^3 + {p^ + ^'^74 



di + A{p^ + 3g2)d3 + 2^3(p2 + 3g2)d2 + i8g(p2 _ ^2)^2^3 + (p^ + 3^2)2^2 
-^{2 [di - Qqd2 + 2(3p2 + g2)d3]3/2 

+2 [di + 3(p + q)d2 + 2(/ - Apq + 7g2)(i. 
+2 [di - 3(p - g)(i2 + 2(/ + 4pg + Tg^)^. 



3/2 



3/2 



j3/2 
i3/2 



+ 



+ 



di + -{p^ + 3g2)d3 - -V27(p2 + 3g2)d2 + 549(^2-^2)^2^3 + (p2 + 3^2)2^2 



f^i + + 3g^)d3 + -V27(p2 + 3g2)d2 + 549(^2-^2)^2^3 + (p2 + 3^2)2^2 



3/2 
3/2 



(29) 



where the first term is the classical potential, the second one comes from one-loop vector dia- 
grams and the fourth and the fifth from one-loop scalar diagrams of H and respectiveljU. 
The effective potential has a symmetry arising from the permutations of the diagonal elements 
of the background (M) and has the following invariance: 



Vi{di,d2,d'i;q,p = ±3g) = ^1(^1,^2,^3; -2g,0). 



(30) 



More generally, the potential is invariant under rotations of 27r/3 in the {p, |g')-plane and in the 
reflections of p: 



p-q 
2 

2 



-p-q 
2 

2 



and 



P 



-P 



(31) 



Thus there is a fundamental region, which determines the potential over the whole plane. We 
choose the fundamental region to be bounded by the two lines p = and p = —3q together with 
the condition p > 0. 

In the fundamental region, there can be four different minima at the critical parameter 
values di, d2 and ^3. The one at the origin (denoted by 1 in fig|5]) is the symmetric minimum. 



®By dropping the last term, i.e., the five last lines and scaling di y, d2 ^ ijs, and ds — > 2x, one obtains the 
effective potential for EQCD in the presence of a finite (imaginary) chemical potential using the notation of |27] 
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Figure 2: 1-loop effective potential in the ((7,p)-plane at the critical point for — 0.01 and d2 — (left 
panel) and ^2 = 0.05 (right panel). Light areas represent the minima of the potential. Solid lines separate 
the three identical sectors which are related by the permutation symmetry of the diagonal elements of the 
background field (M). In the absence of ^2 there is an additional U{1) symmetry making the directions 
marked with dashed lines identical to the p — direction. This symmetry explicitly broken by finite c?2 
as seen on the right panel. 




the minima 2 and 3 are connected by the permutation symmetry and correspond to the same 
physical broken minimum, with Tri?^ < and Tr^^ = 0. The minimum 4 corresponds to a 
phase with TrH'^ > and TrA^ = and is connected continuously to the minimum 2 by a global 
U{1) symmetry ii d2 = 0. If ^2 7^ 0, the U(l) symmetry is lost and the minima 2 and 4 are no 
longer equivalent. If ^2 > the minimum at 2 is favored over 4 and vice versa. 

Setting d2 to zero and expanding in d^ up to order 0(^1) the potential reads (for p = 0): 



Vi{du0,d3;q,0) = ISq^ds 



4df^ ,2n 

(32) 

SO that in the limit — > the potential has two coexisting minima and a first order transition 



27rd3j 3(i3 V ^ 47r2d3 



Qtt ) 



for 



35? 



, ^ .J 0.0759909^, (33) 

This sets the scahng of the critical hne as a function of at smah ^3. Corrections to this, and 
d2 dependence, are obtained by minimizing the real part of Eq. (j29|) numerically. The results are 
shown in Figsll]and[5l The phase transition is accompanied with a discontinuity in q: 



Aq=\ 

'Zbrokon 'Zsymmetric | 



2TTd3' 



(34) 



We see that the transition gets stronger as the coupling ^3 gets smaller justifying a posteriori 
the semiclassical approximation. 
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Figure 4: 1-loop perturbative phase diagram of the soft potential, Vi, as function of di,d2 and ^3. A 
first order critical line separates the two phases. The symmetric phase refers to the phase where with 
d2 = the order parameter vanishes and with (i2 7^ is smaller than in the broken phase. The data points 
represent non-perturbative lattice measurements with ^2 = on a ~ 12^ lattice. The perturbative 
result approaches the lattice data points for small values of ^3 where the transition is very strong. Small 
discrepancy between the perturbative result and lattice data points at small ^3 / 173 is mostly due to finite 
volume effects. 



4.2 Lattice analysis 

The perturbative calculation is valid only for small ^3 and a non-perturbative lattice analysis 
has to be performed to obtain the full phase structure of the model. For the simulations we 
used a hybrid Monte-Carlo algorithm for the scalar fields and Kennedy-Pendleton quasi heat 
bath and full group overrelaxation for the link variables |28t [29} [30j . 
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Figure 5: 1-loop perturbative phase diagram of the soft potential, Vi, as function of di and d2 with 
da = 2. A first order critical line separates the two phases. The non-analyticity ai d2 = is due to the 
change of global minimum between minima 3 and 4. 



The transition was found to be of the first order for ah parameter values used in the simula- 
tions ((^3 < 4 and d2 < 0.15) accompanied with a large latent heat and surface tension; hysteresis 
curves showing discontinuity around critical point in (TrM^M)j^ can be seen in FiglHJ The 
probability distributions of TrM^M along the critical curve are very strongly separated (see 
FigJT]). This makes the system change its phase very infrequently during a simulation, and mul- 
ticanonical algorithm is needed to accommodate a phase flip in reasonable times for any system 
of a modest size |31j . Even with the multicanonical algorithm, the critical slowing restricts us 
to physical volumes np to V <50/g^. 

The pseudo-critical point was determined requiring equal probability weig ht for TrMtM m 
both phases. The simulations were performed with (3 = 12 and a lattice size = 12^, which 
precludes the continuum extrapolation as well as the thermodynamical limit. However, these 
limits were studied for one set of parameter values and the dependence of the critical point on 
both lattice spacing and volume were found to be of order of five per cent for the lattice spacings 
and volumes used (see Figl6]and Tabled]). 

The phase diagram can be seen in Figl9]and FigllOl The non-perturbative critical line follows 
the perturbative one for small values of ds, but for larger fluctuations make the system prefer 
the symmetric phase. The discontinuity in (TrM^M) along the critical line diminishes, as gets 
larger (see Fig Jlip . but it seems that the discontinuity persists, even if its magnitude diminishes 
in the limit d^ ^ oo suggesting that there is a first order phase transition for any (positive) 
value of da. 

5 Conclusions 

In this paper, exact relations between the lattice and continuum MS regulated formulations of 
the Z(3)-symmetric super-renormalizable effective theory of hot QCD, defined by Eqs.([T]),(jl]), 
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Figure 6: Volume dependence of the pseudo-critical point with = 2 and d2 = 0.1. The pseudo- 
critical point was determined by requiring equal probability weight for TrM^M in both phases. The line 
represents a linear fit. The dependence on lattice spacing and volume seem to be within 5% for the lattice 
spacings and volumes used. 



/3 


Lattice volumes 


12 


8^, 10^, 12^ 16^ 


16 


12^, 16^ 20^ 


20 


l6^ 20^, 243 



Table 1: Lattices used in the continuum and thermodynamical extrapolation of the critical point, seen 
in FigEl 



and ([5]), have been calculated. The Lagrangians and the operators up to cubic ones have been 
matched to 0{a^). These results make the non-perturbative lattice study of the theory possible. 

An interesting model with non-trivial dynamics is obtained by setting q = in Eq.(jl|). The 
model amounts to a natural generalization of EQCD to complex variables. The phase diagram 
of the model has been determined using lattice simulations. Two distinct phases were found, a 
symmetric phase with small (TrM^M)j^, and a broken phase with large {TrM^ M)y^. The two 
phases were found to be separated by a strong first order transition with a large surface tension 
and discontinuities in the operators. In contrary to EQCD, where the first order line terminates 
at a tricritical point, the model seems to have a first order transition with all values of ^2 and 

In the future, it is our goal to map out the phase diagram in the full parameter space of the 
model, rather than in a restricted region as in the present exploratory study, in order to search 
for regions in which the phase diagram would resemble that expected for the finite-temperature 
SU(3) pure Yang-Mills theory. 
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Figure 7: Histograms of TrM^Af in logaritlimic scale witli c?2 = along the critical curve. Transition 
channel between the peaks weakens and the transition gets stronger for decreasing ^3. For d^/g^ = 0.5, 
the relative probability density in the tunneling channel is suppressed by a factor ~ 10""'^°. 
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A Details of renormalization 

In this appendix, we give details of the calculation of the renormalization of the mass parameters 
ci and di and the condensates (^g^TrZtZ), (c/g^TrMtM), {g^'^TvZT^Z'^) ,{g^'^2DeiZ) , and 
(53"^2ReTrM3). 

The renormalization calculation compares ultraviolet properties of the two regularizations 
and thus it is irrelevant, in which phase we carry out the computation. We chose to work 
around the symmetric vacuum, since the Feynman rules are the simplest this way. However, 
in this vacuum, all components of the gluon are massless and one therefore has to deal with 
infrared divergences. The infrared divergences in the two regularizations are the same, and 
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Figure 8: Discontinuity in the quadratic condensate in continuum regularization (Tril/^Af )^|g for c?3 = 
0.1, 1, 3. The phase transition gets weaker as the couphng da grows. The metastable regions shrink and 
the discontinuity diminishes. 



cancel exactly in the final results. 

Using the expansion ([8]), the potential 1^ is a function of H and A only. The unit matrix 
commutes with any SU(3) matrix, and the interaction with the gluon field arises from commu- 
tator in the covariant derivative in adjoint representation. Thus the gluons couple, on tree-level, 
only to H and A 

Using this expansion, there are four tree-level mass terms in the Lagrangian that require 
renormalization: 

The coefficients ci and di have to be adjusted such that both regularization schemes give the 
same physical masses for all fields (p, x, H and A. The theory is super-renormalizible, and there 
are divergent contributions up to two-loop level only. The masses are obtained from the low 
momentum properties of the two-point correlator: 

hm ((<A(fc)<A(p))) =<5(=^)(A; + p)-2^ (36) 

The difference between the correlators in the two schemes is in the mass and the wave function 
renormalization. However, the effect of the wave function renormalization is of order 0{a) and 
can be neglected. To get the same masses in the different schemes, we enforce the condition that 
the two-point correlators give identical values in the low momentum limit. The zero-momentum 
lattice correlator can be written in a weak-coupling expansion 

{{mm))a = [{mm)a - {mm)m] + [{mm{-si))a - {mm{-si))m] 
+ i [immsha - immshus] + {{mm))Ms + (^(«)' (37) 

where double and single brackets represent exact and Gaussian expectation values, respectively, 
and subscripts give the regularization scheme. From here, we can read the condition for the 
lattice mass term ci by requiring that the exact correlators give the same value up to order 
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Figure 9: The phase diagram of the soft potential as a function of di,d2 and ^3. First order critical 
line separates two phases. Solid lines represent polynomial fits to the lattice data points and dashed lines 
are the perturbative predictions. The symmetric phase refers to the phase where with d2 ~ the order 
parameter A vanishes and with 7^ is smaller than in the broken phase. Also (TrM^Af) is significantly 
smaller in the symmetric phase and the critical line was determined requiring equal propability weight 
for (TrA/tM) in both phases. 



^, 1 / r 

-Sl))a,lPl - ('/"/'(-5'/))mS,1PI 

^'^'5?)ms,ipi] ) + 0{a 
Similarly, we get for the mass term of adjoint fields: 



- ^1 ^ 1 
ci = — + — 

93 93 
1 

+ 2 



?)a,lPI 





di 1 / 


— - Cl 


+ - + - 


V9l J 


93 91 V 



{HH{-Si))a,rPi - {HH{-Si))^ ypi 



+ 2 



(38) 
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14 



0.5 




-0.2 



-0.1 







Broken 



0.1 



0.2 



Figure 10: The the phase diagram as a function of with d-^ ~ 2. The symmetric phase refers to the 
phase where with d2 — ^ the order parameter vanishes and witli 7^ is smaller than in the broken 
phase. 



The correlators in both regularizations are in two-loop weak-coupling expansion infrared-divergent 
quantities. However, since the infrared properties of the two regularization schemes are the same, 
the infrared divergences cancel exactly in the difference. 

The renormalization of the condensates is done very similarly. The condensates can be 
expressed, in both regularization schemes, as derivatives with respect to mass parameters of the 
free energy and thus they can be related. For the quadratic condensates we get: 



and for the cubic: 



(TrMtM)^ = 

(2ReDetZ)Hs 
(2ReT¥M3)5^ 



MS 



dciiji) 

^fus 
ddi{^j) 

_ ^fus . 
dc2 

ddo 



-- {TTZ^Z)a + 

-- {TiM^M)a + 

(2ReDetZ)a + 
(2ReTrM3)„ + 



djfMS - fa) 

dciifi) 
^(/ms - fa) 

djfus - fa) 
dc2 

OUms - fa) 

ddo 



(40) 
(41) 

(42) 
(43) 
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Figure 11: Discontinuity in (TiM^M) along the critical line d™*((i2,d3)- Dotted lines represent second 
order polynomial fits to the data, and the points on the y-axis represent extrapolations to infinite d^. 
The large c?3 extrapolation yields a finite value suggesting that the transition remains of first order even 
at large d^. 



Due to the super-renormalizability the difference in free energy is dimensionally of the form: 

/mS ~ /a 



1 

47r 



^1,0^ + ^1,0 r ^1,1 



+ 



1 



(47r)2 L 



a 

^2,0—^ + (-^2,1-^ + ^2,2 — 



+ ii^2,0C2 + -^2,1(^2 + -^'2,202^2 

+ -02,053^1 (/^) + -C'2,lC3Ci(/Z) + -D2,2d3Cl(A) 



+ 
+ 



1 



(47r)3 
1 

(4^ 



5 



93 ^r< 53C3 , ^ 5if^3 

3,0 r 03 1 h U3 2 



+ (^3,3— + C34— + C3.5 
a 



4 

a 



^4,093 + C'4,l93C3 + C'4,253'^3 + C'4,353C3 + (^4,393^3 + C'4,4S'3C3C?3 



,2 .2 



C4 5C3 + Ci Qdl + C4 6c|d3 + (74.703^3 



+ 0(a), 



(44) 



where the dimensionless coefficients Aij, Bij, Cij, Dij, and Eij are functions of a dimensionless 
combination afl only. The coefficients Qj- and Dij follow from an z-loop computation. For the 
quadratic and cubic condensates we need to know coefficients Dij and Eij in order to obtain 
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the matching of the condensates to order 0{a^), which follow from a two-loop calculation: 



^1,0 

Di,i 

D2,l 
D2,2 

D2,3 

^2,4 
D2,5 

-E'2,0 
-E'2,1 
E2,2 



-2N,{N- 



-2N,{N- 



, 6 ^ S 
ln—+C+- 
a/x 4 



, 6 ^ S 
In — + C + 7 



1) + 



1 



-12 



1 



In— + C 
an 

InA+c 



In— + C 



(45) 



For the quartic condensates, however, the coefficients Cij are needed and a four-loop lattice 
perturbation theory calculation is required for the matching. For the gluon condensates, also the 
Bij are needed. The coefficients i?2,o and B^^ have been calculated in and respectively. 
The coefficient i?4 has been calculated for A'c = 3 using stochastic perturbation theory in [3l] . 
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B Feynman rules 

Using the expanded fields, the potentials become: 



+^Ac2Tt[H ■H-H] + ^xTr[A ■ H] - ^4c2Tr[A ■ A ■ H] 
+<7i{i2c30^ + l2c3/ + ^^C30V 

+ ^4c302^[i/ • H] + ±^jc3x'Tv[H . H] + ±^^cs^^Tv[A ■ A] 
+ ^4c3x'Tr[^ ■ ^] + ^C3(/.xTr[A • if] 
+^iV6c3(PTr[H ■H-H] + i^4V6c3xTr[A • A • A] 

^M'^^^'^'^'f^ • ^ • ^] + ^^V6c3xTr[A • if • ii] 
+^24c3Tr[A • A • A • yl] + ^24c3Tr[ii' ■ H ■ H ■ H] 

+^16c3Tr[A . A • if • ii] - ^8c3Tr[A • ii • ^ • ii]| (46) 



and 



Vi{Z) =gl [diTv[A ■ A] + diTv[H • H]} 

+gl {-6d2Tr[A • A • ii] + 2d2Tr[ii ■ H ■ H]^ 
+gl^d3Tr[A ■A-A-A]+ AdsTriA -A-H-H] 

-2d3Tr[A -H ■ A-H]+ d-iTi[H ■ H ■ H ■ H]^. (47) 

The gauge part of the scalar Lagrangian in Fourier space (momentum conservation, all 
integrations over Brillouin zone, with measure j^n Jq^i ^'^'^ sums understood) becomes |35j : 

Sz = if {P^q).A'^{p)A\q)A'r{r) + ^.2^^- f^^^^p-j ),5.^A^{p)A%)A'r{r)A'^{s) 

+ if {f^qhH'^{p)H\q)A^Xr) + ||2r-/'''^( )A,H^{p)H\q)A^{r)A^{s), 

(48) 

where we use a compact notation: 

api ~ 2 . api ^ 4 . 2 aPi . 

Pi = cos , Vi = - sm , and p = — > sm 49 

~ 2 ' a 2 ' ^ a? ^ 2 ^ ' 



18 



In addition to these there is the pure gluon and gauge fixing sector [36] 
Sw = lfAU-p)Af{p) +fcf{p)c-{p)+ig^rK_^p,cf{p)c\T)A\{q) 

-^gWir'f'' + r'''fns^P^-C%p)Ai{q)Ai{r)c\s) + ]^gl^^A'i{-p)At{p) + S3 + ^4. 

(50) 

Contributions of three and four gluon vertices 6*3 and ^4 can be found in [37], Eqs. (15. 39), (15. 43) 
and (15.53), where one needs to replace {\){5ab^cd + • • •) with {-^){5ab5cd + • • •) [36j. 



C Calculation of the diagrams 

The perturbation theory calculations were done using symbolic manipulation language FORM [38], 
For formalized computation, it is advantageous to write all the color tensors in the fundamental 
representation, i.e. using the generators of the group: 

Tr [T"T^T^] = ^{(f^^ + ip^^) (51) 

Tr[T"T^] = (52) 

Then all the color contractions in loop calculations can be done systematically with repeated 
use of the Fiertz identity: 

TtjTki = -^{SiiSjk - jj-Sij6ki). (53) 
The following combinations are found in the action: 

•^abc ^ 2Tr(^"[T^^^]) (54) 

^abc ^ 2Tr(^'^{^^ T^}) (55) 

fabcjbde ^ _2Tr[r", ^'=]Tr[^^ T'^] (56) 

^abe^cde ^ 2Tr{r'' , r''}Tr{r^ T"'} - ^(5"^^"' (57) 

In lattice perturbation theory, the numerators of the integrals contain complex trigonometric 
objects. These can be systematically reduced to squares of sines, which also appear in the 
denominator, by repeated use of the following formulae (no summation over repeated indices): 

x + yi= Xiyj + jjiXj (58) 
{x + y)i = xjyj - —xiyi (59) 

x/ = Sii - —xj (60) 
XiyiXjyj = ^({x + y)i - xjxjyj - yfxjyjj . (61) 
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This procedure generalizes trivially also to higher order loop calculations. The set of integrals 
can be further reduced by applying a trigonometric identity (for j > 2) : 

4 J_n (27r)3 (£2 + ^2)i j _ I J_^ (27r)3 {x^ + m2)i-i ^ J_. (27r)3 {x^ + m2)J ' ^ ^ 



D Diagrams for mass renormalization 

In this section, we give the zero momentum diagrams that affect the mass renormalization. 
The expressions are in lattice regularization and the symbol "MS" refers to the result of the 
corresponding diagram in the MS regularization. Solid and wiggly lines represent scalars and 
gluons, respectively. Symbols in parentheses indicate fields running in the internal scalar lines. 
The symmetry factors are included in the coefficients. 

The following diagrams with zero incoming momenta contribute to the renormalization of 
the mass term ci of (f)-&eld (with external lines 0): 

• One-loop: The mass rrii refers to the mass of the field running in the loop. In the difference 
between continuum and lattice regularization, the mass dependence cancels. 



(x) 

(A) 
3 "=^3 47ra 




(63) 



• Two-loop: 

— Terms proportional to c!^gi,: Masses mi, m2 and in the denominator refer to the 
masses of the internal lines and = (ci + di). 



e 



{xx4>) 

{AA<^) 
{HH4>) 
ixAH) 
(AAH) 
I (HHH) 

X H{mi, 1712, ms 



2(iV2 - 1) 

-liV,(iV2 - 1) 
- 1) 



' 9 3 9 ^ 3 



2 at3\ ^3^3 



167r2 



log — + C 



+ C(a) + MS 



(64) 



Terms proportional to csg^: The tadpoles are cancelled by the 1-loop counter terms. 



{A) : -lN,{N^ - 1) 
(H) : -iV,(Ar2-l) 



4' 
53 C3 
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2H{md,md,0) + (2/(0) - /(m^)) (^-^^2/(771^)) + 4m^i/'(md, m^, 0) - a^G{md,md) 



3 ^ V 167r2 



6 . - 
log — + C - <5 



+ 0{a) + tadpoles. + MS 



(65) 



(A) 



iAr,(iV2 _ 1^ 
7V,(7V2 - 1) 



6353 



3/(0)(-5^2)/(mrf) + ^/(O) [-d^.I{md)) - ^mi{md) 



= - ( ^N,{N^ + o{a) + tadpoles. + MS 



(66) 



The following diagrams with zero incoming momenta contribute to the renormalization of 
the mass term di of iJ-field (with external lines H): 

• One-loop diagrams: 



a 



(x) 

(A) 



-C3 

-(2%e-lW(fi3 + 4) ^^3^^"^^) 

-(2A^e-3/Ar,)(c3 + 4) 



--c3-4(Ar,-i/Ar,)(c3 + d3) 



53 



47ra 



+ 0{a) + MS (67) 



= -3glNJ{0) = -S^fA^c:^ + 0(a) + MS 



47ra 



(68) 



=9lNJ{0)=glN,^ + Oia) + MS 



(69) 



Two-loop diagrams: 
— Terms proportional to g^c^: 



e 



' (#i?) 


2ci 






9H3 




(A^(/>) 






iHH<P) 






ixxH) 


9*^3 




ixAH) 






(AAH) 


-2(l-3/iV2- 


|iV2)(c3+4)' 


. (i^i^i?) 


-6(l-3/iV2- 


iiV2)(c3 + J3)^ ^ 



> g^H{mi,m2,m3) 
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(- - 4iV,)ci + {ANl + 24/iV2 - 8)(C3 + 4)' 



167r2 



log — + C 



+ 0{a) +MS 
(70) 



Terms proportional to 93(03 + ^3): 



(A) : -2(iV2-i) 



(^) ■■ -m. 



c 2> 



2/7(md, ma, 0) + (2/(0) - /(m^)) ( -d^,J{md) ) + 4m^i/'(md, m^, 0) - a^G(mrf, m^) 



-8(iV,2-l)(c3 + (i3 



5| 



6 

\og — + C-5 



(A) 



+ 0{a) + tadpoles. + MS 



(71) 



-2(iV2 _ 1) 
-2(iV2 - I) 



93'(c3 + fi1) 



3/(0)(-a„2)/(m,) +a2 <^ ^/(O) (-^^^/(mrf)) - -l(0)/(m,) 



-2(iV2 - 1)5|(C3 + d3 



167r2 



+ 0(0) + MS 



(72) 



Terms proportional to g^: The coupling of the adjoint fields A and H is exactly the 
same as as in EQCD, so the term proportional to can be taken from EQCD |36j . 
However, at two-loop level diagrams with an adjoint scalar loop contribute two times 
since there are two adjoint fields. The diagrams with adjoint loops are the following: 



a 



2yi 



E. (27r)3 (27r)3 x'^x'^{lp + 



a" 



d^x d^y 



iL (27r)3 (27r)3 x'^x'^{y^ + m^) 



+ 



a 



d^x d^y 1 
(27r)3 (27r)3 x^x^ 

d^x d^y 



(27r)3 (27r)3 x'^x'^{]p + m^) 



(73) 



(2x + y)^ 



(X^ + m^)y2y2(-('2. _j_ y-) + 77T,2) 
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N^gl +-H{md,md,0) 



« d^x d^y 1 



i (27r)3 (27r)3 x^x^ jp + 



+ 2m^ 



d^x d^y 1 



1 



1 



.2 z-^- 



(74) 



2 (2^)3 (27r)3 [(~)2p(~2 + ^2)(y2 + ^2) 

The last line can be written in a more familiar form using the definition of p, Eq. (j93l 
and the trigonometric identity Eq. 



a dcx dc'y 



f (2vr)3 (27r)3 [(~pp(~2 + ^2)(y2 + ^2) 

a d^x d^y 1 



4^2a2 (2^)3 (27r)3 J2J2 ^2 _y_^ (27r)3 (27r)3 J2 y2y2 



? d^x d^y 1 1 



+ 0(a) 



(75) 



The infrared divergences in these two diagrams cancel and the sum of the diagrams 
becomes: 



Ar2^4 1 ^1 
^^^167r2 



log — + C 



+ 2p] + 0{a) + MS 
(76) 



There are also two other diagrams: 



E,5?(2!/)', 



(27r)3 (27r)3 x2x2(x2 + m^)(y2 + m^) 



(77) 



2 4I d^x d^y 



J2^ x^{x + 2y)■J2iXj{x + 2y). 



2./ ^ (27r) (27r)' 2,2^2^2;2 -)_ 7722-^^^2 _|_ j^2^|-^^ _|_ + ^d) 



After repeated use of Eqs. ()61|) . the both diagrams can be written in the form 



2 4/ ^ d'x dh x'-2iE^^^yf 



_7L (27r)3 (27r)3 x2x2(x2 + m?^ijp + m^) ' 



(78) 



(79) 



with the negative and the positive sign coming from the first and second diagram, 
respectively, so that their sum cancels exactly. 
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The sum of diagrams proportional to g'^ reads: 



9lK 



167r2 



^2 3A/^,2 



log — + C 
an 



+ 0(a) + MS 



(80) 



It is noteworthy that the scale dependence from the diagrams containing only gauge 
interactions with a single adjoint scalar field cancels exactly in the renormalization. 
However, upon the inclusion of another adjoint scalar field this property is lost. 



E Diagrams for operator renormalization 

In this section, we give the results for the vacuum diagrams that affect the renormalization of 
quadratic condensates present in the action. The diagrams needed for the quadratic condensates 
(TrZtZ)^^, (TrMtM)j^, and {TrZ^TvZ}^: 

• One-loop: 



O 



X 
H 
A 



ij(m) 
ij(m) 

l{N^-l)J{md) 
i(7V2-l)J(m<i) J 



} = (N^ - 1) J(md) + Jim) 



(81) 



Two-loop: the counter terms cancel the {ca, c/aj-dependent linearly divergences terms, and 
only the gauge diagrams contribute: 



oo 



H 
A 



I) 



-93 



^I(O) - 6I{0)I{md) - a^mlliO)I{md] 



(82) 




A 
H 



-gl[-I{md)Iimd) + 4I{0)I{md) 



- imlH{md, rrid, 0) - a^G{md, m^)] 

The diagrams needed for the cubic condensates: (2ReDet.Z) and (2ReTrM^): 
• Two-loop 



(t)hh 
(paa 
Xah 
hhh 
aah 



18^2 

6^-2 
1 fAT2 



c 

t2 



l)cl 



— (^c 

f{l/N,-\N, + \Nl){\c2 + d2? 
%l/N,-lN, + \N^){\c2 + d2f 



(83) 



> g^H{mi,m2,m3) (84) 



+ ^ 1 + 12 ( 1/iVe - ^iVe + (1c2 + ^2) 



glH{m,,m2,ms) (85) 
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F Basic lattice integrals and numerical constants 

In this appendix, we list the basic lattice integrals and numerical constants defined and calculated 
in [Ml Eg [36]. 

Integrals: 



J(m) = 



I{m) 



H{mi, 1712,1713) 



(2^ 



ln(x^ + 



H{0) ^ 1 



d^x 



Sm2 



(27r)3 + m? 



1 



m + 0{am'^ 



d^x d^y 



(27r)3 (27r)3 x'^ + mlP + ml + ^2 



167r2 



^ — r + C + ^ + 

a(mi + m2 + ma) 2 



G{m, m) = 



In 

d^x d^y 



f (27r)3 (27r)3 (~2 + ^2)(y2 + f„2)(;^2^ 



^4 



6 m? 
4^ 



H'{mi,m2,m3) = {d„a)H{mi,m2,m3) 



(86) 



(87) 



(89) 
(90) 
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Numerical constants: 



1 r/^ o 



Ki 



K4 



1 r^n 



27r4 



3.175911535625 



^^sm'^{xi)sm'^{xi + m 



-tt/2 
1 f^/^ 



P = 



n/2 



d^xd^y 



(Ei sin2(xi))2 sin2(yj) sin(xfc + y^) 

,2/^ \ 2 



1.942130(1) 



(91) 
(92) 



X^jSin (xi)sm (xj +yj) sin (x) 



-0.313964(1) 

1 /■V2 



47r4 



d^xd^y 



(Ei sin (xi))2 ^. sm{xj + yj) sin (xi)) 

J2i sin2(xi) sin2(xj + yi) 



(E.sin^(y.))^ 



-7r/2 



1 /■V2 



d^xd^y 



C = lim 

2^0 



7r/2 
1 /•'^/2 



47r4 
0.08849(1 



sin^(xi) sm^ivj) sin(xfc + y^) 

. sin^(xi) sin^(xj + yj) sin^(yt) 
(^^ sin2(xi))2 ^. sin2(yj) sin(xfc + y^) 

d^xd^y 



0.958382(1) 
a 1.204295(1) 



(93) 
(94) 

(95) 



tt/2 



(Ei sin2(xi) + z) sin^iyj) Y^^, s\n\xk + yk) 



In 



3 1 

z ~ 2 



(96) 
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